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ABSTRACT 

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near 
the reionization epoch at z = 6.595. The giant LAE is found from the extensive 1 deg^ Subaru 
narrow-band survey for z = 6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and 
subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 
LAE candidates, this LAE is not only the brightest narrow-band object with L(LyQ;) = 3.9 ±0.2 x 10'*'^ 
erg s~^ in our survey volume of 10^ Mpc'^, but also a spatially extended Lya nebula with the largest 
isophotal area whose major axis is at least ~ 3". This object is more likely to be a large Lya 
nebula with a size of > 17-kpc than to be a strongly-lensed galaxy by a foreground object. Our 
Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is 
wfwhm = 251±21 km s^^, and that the line-center velocity changes by ~ 60 km s~^ in a 10-kpc range. 
The stellar mass and star-formation rate are estimated to be 0.9 — 5.0 x IO^^Mq and > 34 AfQyr"^, 
respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep 
Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this 
could be an important object for studying cooling clouds accreting onto a massive halo, or forming- 
massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment 
of inter-galactic medium. 
Subject headings: galaxies: formation — galaxies: high-redshift — cosmology: observations 



^ Based in part on data collected at Subaru Telescope, which is 
operated by the National Astronomical Observatory of Japan. 

^ Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 

^ Observatories of the Carnegie Institution of Washington, 813 
Santa Barbara St., Pasadena, CA 91101 

^ Carnegie Fellow; ouchi _at_ ociw.edu 

^ Department of Astronomy, School of Science, University of 
Tokyo, Tokyo 113-0033, Japan 

'' Department of Astronomy, University of Arizona, 933 N Cherry 
Avenue, Rm. N204, Tucson, AZ 85721-0065 

Research Center for Space and Cosmic Evolution, Ehime Uni- 
versity, 2-5 Bunkyo-cho, Matsuyama 790-8577, Japan 

* Kavli Institute for Particle Astrophysics and Cosmology, Stan- 
ford University, 2575 Sand Hill Road, Menlo Park, CA 94025 

^ Department of Astronomy, Cornell University, Ithaca, NY 
14853 

^^ Astronomy Centre, University of Sussex, Falmer, Brighton, 
UK 

^^ Optical and Infrared Astronomy Division, National Astro- 
nomical Observatory, Mitaka, Tokyo 181-8588, Japan 

^^ Nobeyama Radio Observatory, Minamimaki, Minamisaku, 
Nagano 384-1305, Japan 

^^ Subaru Telescope, National Astronomical Observatory of 
Japan, 650 North A'ohoku Place, Hilo, HI 96720 

Astronomical Institute, Graduate School of Science, Tohoku 
University, Aramaki, Aoba, Sendai 980-8578, Japan 

^^ Department of Physics and Astronomy, University of British 
Columbia, 6224 Agricultural Road, Vancouver V6T IZl, Canada 

^^ SUPA Institute for Astronomy, University of Edinburgh, 
Royal Observatory, Edinburgh EH9 3HJ, UK 

Centre for Astrophysics, Science & Technology Research In- 
stitute, University of Hertfordshire, Hatfield ALIO 9AB, UK 



1. INTRODUCTION 

Identifying the first stage of galaxy formation is one of 
the ultimate goals in astronomy today. Theoretical mod- 
els predict that primordial gas accretes onto the center 
of halos via gravit ational cooling with subsequent star - 
formation activity (jFardal et al.ll2001l : lYang et al.ll2006D . 
These primordial galaxies make spatially extended Lya 
nebulae caused by hydrogen cooling, and it is suggested 
that high-z extended Lya sources, o r Lya blobs, are can- 
didates for p rimord ial galaxies (e . g. iMatsuda et al]|200^ 
Saito et al.l l2006t iNilsson et atl l2006t ISmith fc Jarvii 
2007t l. Lya blobs are found mostly at z ~ 2 — 3, 
and have angular extents of ~ 5 — 16 arcsec with to- 
tal Lya l uminosities ranging from ~ 6 x 10''^ to 10^'' 
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The most prominent 



Ly g nebulae kii o wn to date are blobs 1 and 2 found 
bv ISteidel eFai] (|2000D . which extend over > 100 kpc 
with L{Lya) ~ 10''^'' erg s~^. Although Lya blobs are 
candidates for galaxies with gas inflow of cooling accre- 
tion, it is also suggested that Lya blobs can be pro- 
duced by inten sive starbursts ass ociated with significant 
outflo ws (e.g., iTaniguchi fc Shio va 2000; Wilman et aLJ 
I2005f l. by a hidden AGN (e,£^Haiman &Reeai200li), 
or by both of them (e.g., iDev e t al. 200l; Yang et al. 
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Fig. 1. — Optical to infrared images of Himiko. North is up 
and east is to the left. We display 10" X 10" images at BVRi' z' 
and NB {NB921) bands from Subaru/SXDS, at JHK bands from 
UKIDSS-UDS DR3, and 3.6 - 24/.tm bands from Spitzer/SpUDS. 
We show intensity contours in Af_B921 and z' images. The black 
contours denote 1 a level of sky fluctuation. The yellow contours 
represent (2, 3, 5, 10, 15) and (2,3,4) a levels of sky fluctuations 
in NB921 and z' images, respectively. We also plot the position of 
the DEIMOS slit by the red box. The dispersion direction towards 
red spectrum is shown by the magenta arrow. The green arrows 
point to the position 1 and 2 that are probable peaks in the z' 




Fig. 2. — Composite pseudo-color image of Himiko. The RGB 
colors are assigned to 3.6/im , z' , and NB921 images, respectively. 
North is up and east is to the left. The image size is 5" X 5". The 
white bar at the bottom right represents the length of one arcsec- 
ond. The brightest peak with a bluish white color corresponds to 
position 1. The position 2 is located 1.1 arcsec west of the position 
1. 



in preparation). In fact, the i nfrared-subm ni and X- 
ray observations indicate that ISteidel et aLl 's blobs 1 
and 2 would be powered by a heavily obscu red star- 
burst (iGeach et al.ll2007t iMatsuda et all l2007l: see also 
Chapman et al.ll2004D"and an AGN (iBasu-Zych fc Scharl 
2004() . respectively. IMatsuda et all ()2006h claim that 
all of their spectroscopically-identified Lya blobs are 
likely to be the sites of massive galaxy formation be- 
cause of their large line widths of wfhwm ^ 500 
km s~^. It is also well known that such bright 
large Lya nebulae are associated with radio-loud 
AGN (e.g. iMcCarthv et al.lll987l: Ivan Ojik et "all 119971: 
Reuland et all 120031 : iBarrio et al" '2008': 'Smith et all 
20091) or radio-quie t quasars (Wcidingcr ct al. 200J!^; 
Hennawi et al.l l2008f ) . In this way, extended Lya neb- 



ulae shed light not only on primordial galaxies but also 
on massive-galaxy formation and AGN activity. 

So far, extended Lya nebulae have been fo und only at 
z = 2 — 5 with the majority at 2; ~ 2 — 3 (e.g. , I Saito et all 
|2008( ). Due to this current observational limit, it is diffi- 
cult to identify primordial galaxies as well as to study the 
early stage of massive galaxy formation. At z ~ 2 — 3, 
the mean metallicity of the inter-stellar mediu m is al- 
ready as high as Z = O.IZq (jSadat et al.ll200ll : see also 
iPettini et al.1 [l997r ). It is predicted that the fraction 
of primordial galaxies to met a.1 enriched galaxies woul d 
be quite low at z ~ 2 — 3 (|Scannapieco et aP l2003f ). 
and that the fraction of population III to population II 
star-formation rat e (SFR) rises with increasing redshift 
(jTrac fc CenI 120071 ). On the other hand, the importance 
of the early stage of massive-galaxy formation has been 
recognized by the downsizing behavior of stellar-mass as- 



sembly (jCowie et al.lll996l ) . A massive population at z = 
2—3 selected from a distant-red galaxy sample is old, 2—3 
Gyr, and their mean format ion reds hift is estimated t o 
be z > 5 (|Labbe et al.l 120051 : see also lKriek et"ani2006[ ). 
It is implied that galaxies at the massive end would have 
a very high specific star -formation rate (SSFR) at z > 4 
(jDrorv fc Alvareal2008l ). and that a major active star- 
formation in massive galaxies probably takes place at 
z > 4. Thus, it is important to study extended Lya neb- 
ulae at a redshift greater than the current-observational 
limi t, especially at the reionization epoch of z ~ 6 — 
11 (|Fan et al.|[200l iKomatsu et al.l [2008). This epoch 
is a lso the today's observational limi t of galaxy stud- 
ies (live et all 120061: IStark et al.l l2007l : lOta et al.l 120081 : 
iBouwens et al.ll2008[ ) . Moreover, such bright Lya sources 
can be a good laboratory for understa nding reioniza- 
tion (e.g. [K ashikawa et al. 2006; Diikst ra et all l2007bl : 
iMcQuinn et al., 2007; Kobavashi ct al. 2( J07|) and meta l 
enrichment (e.g. Martin ct al. 2002: Bouchj et al.ll2007l) 
of inter-galactic medium (IGM). 

In this paper, we report our discovery of an extended 
Lya nebula, which we named Himiko ^^ , near the reion- 
ization epoch at z = 6.595. We describe the photo- 
metric identification and spectroscopic confirmation of 
this object in fj2l and present detailed properties such 
as kinematics and stellar population in ^31 We discuss 
the nature of this object and prospects of future ob- 
servations in f|4l Throughout this paper, magnitudes 
are in the AB system. We adopt {h,flrm^A,ns,<Ts) — 
(0.7,0.3,0.7,1.0,0.8). 

^^ Himiko is a name of legendary queen in ancient Japan. 
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Fig. 3. — Surface brightness (SB) profiles of our object. Open 
squares show a SB profile in A'^iJ921 image, and filled circles denote 
that of a continuum-subtracted (Lya) image. The dotted line is 
the Gaussian profile representing the point-spread function of our 
image with an FWHM = 0".8. 

2. DISCOVERY 

2.1. Photometric Identification 

We have identified a candidate very bright spatiaUy- 
extended Lya emitter (LAE) at z ~ 6.6 in the 
course of our deep and wide-field narrow-band imag- 
ing program in th e Subaru/XMM-N ewton Deep Sur- 
vey (SXDS) field (jOuchi etal.l 120081 ). In 2005-2007, 
we took narrow-band images in the NB921 filter with 
a centra l wavelength of Ac — 9196A and a FWHM 
of I32 A (iHavashino et al.l[2003[ ) using Subaru/Suprime- 
Cam (jMivazakiet al.ll2002f l. The 1 deg^ field is cov- 
ered by 5 pointings of Suprime-Cam with a total on- 
source integration of 4 5.1 hours. T h ese data are reduce d 
with SDFRED (.Yagi ^eTall [20f)l lOiichi et all [200I . 
and aligned with o ptical broad-band images of SXDS 
(jFurusawa et al.|[200S) . The FWHM of the seeing size in 
the aligned images is ~ 0".8. The 3a limiting magnitude 
in NB921 is 26.2 — 26.4 mag in a 2".0-diameter aperture. 
Combining the deep optical broad-band images of SXDS, 
we have selected candidates of z ~ 6.56 ±0.05 LAEs that 
satisfy our photometric criteria of the narrow-band ex- 
cess (z' — NB921 > 1.0), no detection of blue contin- 
uum flux {B > B2cr and V > ¥20-), and the existence 
of Gunn-Peterson trough {[z' < z'^^ & i' — z' > 1.3] or 
\z' > z'^„]). The B2a and ¥2^ are defined as 2cr limiting 
magnitudes of B and V bands, respectively {B2cr — 28.7 
and V2(T — 28.2), while z'^^ is the 3a detection limit 
(zg^ = 26.5). We have obtained a photometric sample 
of 207 LAEs at z ~ 6.6 down to NB921 = 26.0 in a 
comoving survey volume of 8 x 10^ Mpc^. The sky dis- 
tribution of our LAEs show a rectangular area (8' x 20') 
with a number density of LAEs higher than the aver- 
age by a factor of 2. In this high-density region, we 
find the object, Himiko, that has the brightest NB921 
magnitude and the largest isophotal area among the 207 
LAE candidates. We present snapshot images in Fig- 
ure [1] and a close-up color composite image in Figure 



1 — 1 


b 


1 1 1 


' 




1 1 1 


d 




Himiko 


Himiko 








,^ 


■ 


■ 


-1 - 















Sh 












rrt 














4 


- 




- 


- 


cfl 




fe: 


m ^ 






F-i 


3 


^U^ ^^^ 


JH 






< 


^fH 


et* 


















hH 




'1 ' My, 


XiV^l ¥1 




1 


^ 


? 


'ii^t^rjBrfi Ij ,- 


W^^^i 


^- 


[ 







1 — .M wMft||4^^fiWi 


noiffii^WMB 


rifM 




xi 




'lijNfTffilMl'jhi 


^iJlMlllMlll 


ip 




a 




^^all-f )jjnfffyil|i 


^ipjSHilililff 


&T 









'^^^^g 


.^^|4IUI|yH|iK 


^ 


1 — ■ 


(/J 


— ^1 1^^ ^^^^ 




n 


1 , 1 °v^^ 






\ 




; ° 



23 24 25 

7V592/ (total) 



25.5 26.5 20 

(SB) N 



Fig. 4. — Isophotal area of our 2 = 6.6 LAE candidates as a 
function of total A'^_B921 magnitude (left panel) and average NB921 
surface brightness, {SB) (middle panel). (SB) is expressed in units 
of mag arcsec"^. The square with a label of Himiko is our giant 
LAE. The open squares with a filled square represent 2 = 6.6 LAE 
candidates showing a possibly extended profile with a FWHM of > 
1".2 in the AfB921 image (a PSF FWHM of 0".8), while the simple 
open squares indicate the other (FWHM~ 0".8 — 1".2) candidates. 
The measurements of FWHM include large uncertainties for faint 
sources with 7V_B921 ~ 25 — 26, and become unreliable in this 
faint magnitude regime. The gray star and cross marks denote our 
simulated blob 1 and 28, respectively (see i|4.1l l. The right panel 
plots the isophotal-area distribution of the LAE candidates with 
the thick line. The thin line presents A' = 0. 

H The total magnitude of Himiko is NB921 = 23.55, 
which is brighter than the second brightest candidate 
(iVB921 = 24.06) by 0.5 magnitude. This object is sig- 
nificantly extended, in contrast to the compact point-like 
profiles of the other LAEs. In Figure [3l we plot sur- 
face brightness profiles of this object in NB921 image 
and a continuum-subtracted Lya image obtained from 
iVi?921 and z' d ata with the assumptio n of z = 6.6 LAE 
spectrum shape (jShimasaku et al.ll2006f ). Figure [3] shows 
that the profiles of this object is clearly more extended 
than that of point-spread function (PSF) of our image 
{FWHM = 0".8). The outskirts of the profile reach a 
radius up to ~ 1".5 and a possibly to ^ 2" . If we define 
the isophotal area, Ajgo, as pixels with values above the 
2a sky fluctuation (26.8 mag arcsec"^ in iVi?921), the 
isophotal area of Himiko is Ajgo = 5.22 arcsec^ in the 
iVi?921 image. Figure [3] presents the isophotal area of 
our z — 6.6 LAE candidates as a function of total 7VB921 
magnitude and average NB92\ surface brightness. The 
average A^i3921 surface brightness, {SB), is the value of 
an isophotal flux divided by the isophotal area, where 
the isophotal flux is the one summed over the isopho- 
tal area. We mark a possibly extended (FWHM> 1".2) 
sources with a filled squares to distinguish between bright 
point-like and faint extended sources with a compara- 
ble isophotal area. The right panel of Figure |4] plots an 
isophotal-area distribution with a histogram, which vi- 
sualizes the extended nature of our object among all of 
the LAE candidates. Figure |4] indicates that there are 
no LAEs similar to Himiko. We confirm that the bright- 
est source from the previous 0.2 deg^ Subaru Deep Field 
(SDF) survey for z ~ 6.5 LAEs is only as bright as our 
secon d brightest candidate with no significant spati al ex- 
tent (|Taniguchi et al.ll2005t iKashika wa eral]|2006f) . and 
that our object is distinguished from all the other z ^ 6.5 
LAEs found in the previous studies. By the comparisons 
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with z ~ 6 — 7 LAEs found in previous studies, we no- 
tice that the size of this object is outstanding. The ma- 
jor axes of the isophotal area in NB921 and z' bands 
are ci 3".l and ~ 2".0, respectively. Additionally, the 
NB921 (z') image shows potential diffuse components 
which continuously extend by ~ 1" (~ 0".3) around the 
isophotal area with a surface brightness above la sky 
fluctuation (Figures [1] and [2]). Thus, the size of our ob- 
ject is probably > 3".l and > 2".0 in NB921 and z' 
bands, respectively. Given the fact that this LAE has 
the unusual brightness and size, we refer to this object 
as the giant LAE. 

Interestingly, this object is detected at the 4(7 level in 
the medium deep 3.6/xm image from the Spitzer legacy 
survey of the Ultra Deep Survey field (SpUDS; PL. J. 
Dunlop; Figure [T]), while we find only marginal de- 
tections (~ 2 — 3ct) ^^ in the near-infrared (NIR) im- 
ages from the UKIRT Infrared Dee p Sky Survey T hird 
Data Release (UKIDSS-DR3: Lawre nce et"alll2007( ). We 
ahgn Spitzer/SpUDS and UKIDSS-DR3 images with the 
SXDS optical images, referring a number of stellar ob- 
jects in the field. The relative astrometric errors are 
estimated to be ci^ 0".04, ~ 0".ll, and ~ 0".35 in 
rms, for optical-NIR, Spitzer/IRAC(3.6 — 8.0/xm), and 
MIPS(24/im) images, respectively. We summarize total 
magnitudes/fluxes and 2"-diameter aperture magnitudes 
of Himiko in Table[T] We defi ne the total magnitude w ith 
MAG.AUTO of SExtractor (jBertin fc Arnoutsl[l996l ) in 
the optical and NIR bands. The total magnitudes of 
Spitzer/IRAC and MIPS bands are obtained from a 3"- 
diameter a pertur e and an aperture correction given in 
lYan et all ()2005f ) and the MIPS web page ^3^ respec- 
tively. Note that our object is detected in the 3.6/Ltm 
band, but not in the 4.5/im band. This is probably due 
to the higher noise level in the 4.5/im band, as we ex- 
pect the object to have a fairly flat spectrum at these 
wavelengths. Our measured 3.6/im magnitude of 24.02 
(4cr) would result in a < 3cr detection at 4.5/Ltm for a 
flat spectrum (constant AB magnitude), consistent with 
our tentative 1 — 2cr detection and the large error with 
24.62 ±0.73 24. 

2.2. Spectroscopic Confirmation 

We carried out spectroscopic follow-up observa- 
tions with Keck/DEIMOS and Magellan/IMACS. The 
DEIMOS observations were conducted with the 830G 
grating and the GG495 fllter in the non-photometric 
night on 2007 November 5. The IMACS observations 
were made in the nod-and-shuffle mode with the 150 
1 mm~^ grism and the GG455 filter under the photo- 
metric night on 2007 November 11. The on-source ex- 
posure times of our DEIMOS and IMACS observations 
were 10800 and 15600 seconds in the 0".8 - 1".0 and 
0".5 — 0".8 seeing conditions, respectively. We chose 1".0 
for the slit width in these observations. The DEIMOS 
slit position is plotted in Figure [1] The slit position of 
IMACS is the same as that of DEIMOS, but the posi- 
tion angle is slightly rotated by -1.1 deg. The spectral 

^^ We estimate the 2(t limits of total magnitudes in the vicinity 
of this object to be J = 24.3, H = 24.0, and K = 23 .8. 
^^ http://ssc.spitzer.caltech.edu/niips/apercorr/ 1 
^* See ^3.1l for a possible inclusion of emission lines in the IRAC 
bands. 





TABLE 1 




Photometry of Himiko 


Band 


Mag(2") 


Mag/Flux{Total) 




(1) 


(2) 


/(0.5 - 2kcV)2 




<6 X lO-ii^ 


B 


> 28.7 


>27.9 


V 


> 28.2 


>27.4 


R 


> 28.1 


>27.3 


i' 


> 28.0 


>27.2 


2'3 


25.86 ±0.20 


25.45 ± 0.27 


7VS9213 


23.91 ±0.04 


23.55 ±0.05 


mo. 95 


25.74 ± 0.64 


25.18 ±0.73 


J^ 


24.95 ± 0.53 


24.01 ± 0.43 


i?4 


> 24.7 


>24.0 


K4 


24.42 ± 0.50 


>23.8 


m(3.6Atm) 






24.02 ± 0.27 


m(4.5/^m) 






>23.9 


m(5.8/^m) 






>22.0 


m(8.0^m) 






>21.8 


m(24^m) 






> 19.8 


S{m0^lm) 






< 12mJy 


/(1.4GHz) 






< 100/iJy 



Note. — Col.(l): The 2"-diaincter aperture mag- 
nitude. Col. (2): The total magnitude or flux. In these 
two columns, the upper limits are 2a and 3(7 magni- 
tudes in B — K and 3.6 — 24/iin bands, respectively. 

In units of erg cm~ s~ . 
^ Isophotal magnitudes are 25.67±0.21 and 23.66±0.04 
in z and N B921 bands, respectively. 

Magnitudes in J and K bands are slightly over the 
2a level. However, neither of them are detected beyond 
the 3ct level. 

coverages of DEIMOS and IMACS are 5700 - 9500A and 
4500 — 97OOA, respectively. The spectral resolution of the 
DEIMOS data at 9200A is i? ~ 3600, while the one of 
the IMACS data is i? ~ 700. We have reduced our spec- 
tra with spec2d 25 and COSMOS pipehnes for DEIMOS 
and IMACS data, respectively. Both spectra have a 
strong single line with no detectable continuum. The 
line-center wavelength of the single line is 9232. 7A from 
our DEIMOS spectrum, which coincides with the mea- 
surement from our IMACS spectrum (9233. 5A) within 
lA. Figure [5] presents our spectra in the wavelength 
around this single line. We have conflrmed that the spec- 
tra show no signatures of an [Oiii] 5007 emission line (at 
~ 7044A) and an [On] 3727 emission hne (at ~ 5243A) 
from a z = 0.407 Ha emitter or an [On] emission line 
(at ~ 6872A) from a z = 0.844 [Oiii] emitter, and found 
that this object is neither a foreground Ha nor [Oiii] 
emitter. We cannot distinguish between an [On] emit- 
ter at z = 1.477 and a Lya emitter at z = 6.595 from 
a detection of the other emission line, because our spec- 
tra do not cover a wavelength that would have another 
strong emission line such as Lya, [On] , [Oin] , and Ha. 
However, the DEIMOS spectrum has a FWHM spectral 
resolution of 2. 6 A that would have enabled us to iden- 
tify an [On] AA3726,3729 doublet at z = 1.48 with a 
separation of 6.9A (vertical arrows in the middle panel 
of Figure [5]). Our DEIMOS spectrum confirms no such 
signature of [On] doublet, but a clear asymmetric line 
profile with an extended red wing that is typical for 
a high-z Lya line. We measure the skewness, S, and 

^^ The analysis pipeline used to reduce the DEIMOS data was 
developed at UC Berkeley with support from NSF grant AST- 
0071048. 
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Fig. 5. — Spectra of the giant LAE, Himiko. The top panel shows 
the two-dimensional spectrum obtained from DEIMOS observa- 
tions. The horizontal arrows point to the position 1 and 2. The 
middle and bottom panels present the spectra taken with DEIMOS 
and IMACS, respectively. We show spectra of the giant LAE (solid 
line) and the background sky (dotted line). The units of the ver- 



tical axis are in 10 



in the bottom panel 



(IMACS), and arbitrary in the middle panel (DEIMOS). Two ver- 
tical arrows in the middle panel indicate the wavelengths of [On] 
AA3726,3729 doublet from a 2 = 1.5 [Olll emitter. 



the weighted skewness, S^, defined bv lKashikawa etldl 
11006). We obtain 5* = 0.685 ±0.007 and S'™ = 13.2±0.1 
for our line. Since the average values of z ^ 6.5 LAEs are 
S = -542 ± 0.007 and Su, = 11.5 ± 0.2 (jKashikawa et all 
l2006( l. the line shape of our object is similar to (or more 
positively skewed than) the average. If this line were an 
[On] doublet, the line shape would be negatively skewed. 
Thus, we conclude that this object is a real LAE at 
z = 6.595 with a clear red wing in the asymmetric Lya 
line. This is the spectroscopic confirmation of the giant 
LAE, Himiko, aX z = 6.595. The number density corre- 
sponding to this giant LAE is only 1.2 x 10~^ comoving 
Mpc^'^ aX z — 6.6 We have identified the rare object near 
the reionization epoch. 

We estimate Lya flux, /(Lya), and rest-frame equiva- 
lent width, EWo, to be /(Lya) = 7.9 ± 0.5 x lO'^^ erg 
s-i cni-2 and EWq = lOOtl"^ ^ with z' and NB%2\- 
band photometry in the same manner as lOuchi et al.l 
(|2008f l. The corresponding Lya luminosity, i(Lya), is 
L(Lya) = 3.9 ± 0.2 x 10^^ erg s"^ To check our estima- 
tion, we derive /(Lya) from the IMACS spectrum that 
were taken under the photometric condition. Applying a 
slit-loss correction, we obtain /(Lya) — 11.2±3.6 x 10~^^ 
erg s~^ cm^^ corresponding to L(Lya) — 5.6± 1.8 x 10^^ 
erg s~^. Although the line-flux value from the IMACS 





TABLE 2 




Properties 


OF Himiko 


Quantity 






Measurement 


Redshift(2) 






6.595 


Skewness of Lya (5) 






0.685 ±0.007 


Weighted Skewness of Lya 


(5») 




13.2 ±0.1 


Isophotal Areai(ArB921) 






5.22 arcsec-^ 


Isophotal Area-'^(z') 






1.88 arcsec-^ 


Major Axis2(ArB921) 






3.1" 


Major Kyis?(z') 






2.0" 


LyQ Surface Brightness 




1.51 


x 10~^^erg s~^cm~-^arcsec~^ 


iVB921 Surface Brightness 






25.5 mag arcsec"^ 


/(Lya)^ 




7.9 ± 0.5 X 10-1^ erg s"! cm-^ 


L(Lya)3 






3.9 ±0.2 X 10^3 ergs-i 


M(1250)-* 






-21.67 ±0.73 mag 


Ms 4 






-22.83 ± 0.27 mag 


Lya Line Width (FWHM) 






251 ±21 km s-l 


Rest-Frame EW {EWq) 






lootif A 


Stellar Mass 






3.5t^| X IOIOM0 


SFR from SED fit^ 






> 34MQyr-l 


SFR from UV^ 






25l?^M0 yr-i 


SFR from Lya^ 






36 ± 2Mq yr-i 


Specific SFR 






> 1.6 X 10-^yr-i 


Number Density^ 






1.2 X 10-*^ Mpc-3 



The isophotal areas are defined as pixels with values above the la sky 
fluctuation; 26.8 and 27.3 mag arcsec" in the N B921 and z' images, re- 
spectively. 

The maximum size of the 2ct isophotal area. For isophotal areas above the 
3fT sky fluctuation, we obtain 2' .7 and 1 .2 in the N B921 and z images, 
respectively. 

The Lya flux and luminosity from the photometric measurements. 
* The rest-frame 1250A, M(1250), and B-band, Mb, magnitudes. Since the 
3.6/.im band observes the rest-frame 4180-5168A which is very close to the 
bandpass of B band, no k-correction is applied to the 3.6/im-band magni- 
tude. 

The SFRs estimated from the SED fitting, the UV continuum, and the Lya 
luminosity. The SFRs of UV and Lya are not corrected for dust extinction. 
See i|3.4l for more details. 

The comoving number density corresponding to Himiko. 

spectrum includes the large error, the line flux from spec- 
troscopy agrees with the one from photometry within the 
la error. Since the IMACS spectrum shows no contin- 
uum above the detection limit, EWq cannot be derived 
from the spectrum. We summarize spectroscopic as well 
as various properties of this object in Table [21 

3. A GIANT LAE AT z = 6.595 

3.1. AGN Activity 

We investigate the AGN activity in our object. Our 
spectra show no Nv 1240 line at 9418A as well as no 
broadening of Lya. Moreover, there is no counter- 
part in our Spitzcr MIPS image as well as t he XMM- 
Newton X-r ay and VLA radio catalogs (,Ueda et al.l 
l2008t ISimpson et al.l l2006f ) whose detection limits are 
m(24/im) = 19.8, /(0.5 - 2keV) = 6 x lO'^^ erg cm'^ 
s~^, and /(1.4GIIz) = 100/xJy, respectively. Our ob- 
ject is not also detected in th e SCUBA H alf-Degree Ex- 
tragalactic Survey (SHADES; iMortier et' al. 2005) data. 
We find no significant signal at the location of our ob- 
ject in the submm 850/im map, and obtain the 3cr upper 
limit flux of S'(850/im) < 12 mJy. We present these 
upper limits i n Figure El togethe r with the AGN tem- 
plates used in lOuchi et al.l (120081 ). The AG N templates 
are t he average radio quiet/loud quasa rs (I Elvis et alJ 



119941) ■ and a UV-optical s pectrum of Telfer et al.l (|2002f ) 
and iRichards et aLJ (|2003D that is normalized to the av- 
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Fig. 6. — The upper limits to the spectral energy distribution 
(SED) of our object in the rest frame. The thick bars with an ar- 
row indicate the upper limits that correspond to the VLA (Radio), 
SCUBA (Submm), Spitzer/MIPS (MIR), and XMM-Newton (X- 
ray) results from left to right. Solid and dotted curves are the tem- 
plates of radio-quiet and r adio-loud AGNs taken fromlEl vis et all 
}T99l) . rMfer et al.l ||200^ '). and [Richards ct al. (2003). The thick 
and thin curves represent the template SEDs normalized with hya 
and rest-frame optical luminosities of our object, respectively. 



erage spectral energy distribution (SED) of lElvis et al.1 
We choose the amplitudes (log i/L^) of these 



AGN templates that match to Lya (thick lines) and op- 
tical luminosities (thin lines) of our object. Here, the 
Lya luminosity of the templates is defined by the Lya 
flux that enters into the bandpass of our narrow-band 
(Arest = 1207— 1224A). The optical luminosity is the one 
at rest frame ~ 4700A that corresponds to the Spitzer 
3.6^m band. Figure [S] indicates that our object is prob- 
ably not a typical radio-loud AGN, although the upper 
limit of radio band is marginally below the templates. On 
the other hand, none of our multi-wavelength data can 
place constraints on the existence of radio-quiet AGN. 
Thus, we cannot discard the possibility that our object 
is an AGN. It is also true that there are no positive sig- 
natures of AGN activities in our spectroscopic and multi- 
wavelength data. 

Although it is not clear whether this object has an 
AGN, we investigate how much our broad-band photom- 
etry is impacted by prominent nuclear emission lines in 
the case of AGN. Strong [Oiii] and Ha of AGN would en- 
ter into the 3.6/im and 4.5/im bands, respectively, which 
may boost the magnitudes in these broad bands. We 
estimate the possible contributions from these strong 
fines, assuming the fiux ratios of /LyQ//[oiii] = 4 for 
a type II AGN (|McCarthvlll993l : due to no Lya broad- 
ening ) and fhva/f'Ha = 8.7 for the case B recombina- 
tion |Brockiehurstlll971h . We, thus, obtain 26.3 and 
27.0 magnitudes in 3.6/zm and 4.5/um bands, respectively. 
Since the magnitude of our object is 24.0 in the 3.6/im 
band, the flux contributed from strong lines of AGN 
would be about one-order of magnitude smaller than the 
brightness of our object. Even with an AGN, the mag- 
nitudes of 3.6/ini and 4.5/im bands would include a neg- 



ligible contribution from strong emission lines. It should 
be noted that these flux contributions from strong lines 
could be underestimated, in the case where our Lya flux 
is very strongly absorbed. 

3.2. Possibility of Gravitational Lensing 

Strong lensing is more common among very high- 
redshift objects because of increas ing lensing optical 
depth with increasi ng redshift (e.g.. iHilbert et alll2007t 
iHilbert et al.l [2008) and larger effect of magnification 
bias for very high-redshift objects (e.g.. lWvithe fc Loebl 
120021 ). Moreover, the extended nature of Himiko may 
make the lensing interpretation of this object plau- 
sible. In this section, wc investigate the possibility 
of gravita tional lens i ng. First, we use the catalog of 
Ivan Breukelen et al.l ()2006t ) which shows cluster candi- 
dates at z = 0.5 — 1.5 in this field. The estimated masses 
of these clusters range between 5 x 10"'^^— 3 xIO-'^^Mq. Our 
object is separated from the center of the closest cluster 
candidate by ~ 5'. 8. Due to this large separation, the 
magnification by these clusters is negligible. 

Next, we investigate the possibility of galaxy-galaxy 
lensing. Figure [1] indicates that this source may have 
two peaks in the z' image, and a smooth profile in the 
NB921 image. From a visual inspection, the 3.6/im-band 
profile would appear to be slightly elongated. However, 
our object is not bright enough in the 3.6/im band to dis- 
tinguish between profiles of a real extended source and a 
point source with outskirts made by peaks of background 
fiuctuations. We carry out profile fitting to our z' and 
NB921 images with GALFIT (P eng et al . 200^. We fit 
two profiles whose flux amplitudes and positions are free 
parameters. We find that fitting with two PSF profiles 
leaves large residuals, and that the profile of our object 
is well fit by two circular exponential disks with a half- 
light radius of i?hi = 0".3 which are separated by I'M. 
The positions of two components are determined in the 
z'-band image that shows the possible two peaks, and 
presented with green arrows in Figure [1] On the other 
hand, the two-exponential disk models can reproduce not 
only the z'-band, but also the Af_B921-band profiles. The 
positions of the best-fit models in the NB921 image are 
not the same as those in the z' image, but the differ- 
ences of the positions are only 0".2 — 0".3. Because the 
photometric uncertainties in the z' band are large, it is 
not clear whether two peaks really exist or whether the 
positions of two peaks are different between the z' and 
NB921 images. Nevertheless, we refer to the positions 
of east and west components determined in the z' im- 
age as position 1 and 2, respectively. The brightness 
ratio. Am, oi position 1 (mi) and 2 (7712) components 
are Am = m2 - mi = 0.38 ± 0.38 and 0.65 ± 0.10 in 
z' and NB921 bands, respectively. There are no signif- 
icant differences between Am of z' and iVi?921 bands. 
We cannot reject the possibility that these two compo- 
nents have the same color, and that these components 
are an identical lensed object. Thus, our profile fitting 
does not constrain the possibility of lensing due to the 
large photometric uncertainties of Am. 

Figure [7| presents DEIMOS spectra at position 1 and 2 
whose extraction- aperture size is 0".6 along the slit. The 
skewness and the weighted skewness of the emission line 
are estimated to be {S, 5^) = (0.846 ± 0.018, 17.3 ± 0.4) 
and (0.502±0.023, 9.3±0.4) in the spectra at [position 1, 
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Fig. 7. — Top: DEIMOS spectra at position 1 (blue) and 2 (red). 
The dotted line represents the background sky. The instrumental 
FWHM is indicated with the bar below the emission line. Bottom: 
The line-center offset (upper panel) and FWHM line width (lower 
panel) as a function of position. The AAc is defined by AAc = 
Ac — Ac (Ad = 0), where Ac is the line-center wavelength. The 
definition of Ad is the same as that of the top panel in Figure 
\5\ and Ad = corresponds to the Lya source center. The gray 
open squares are the direct measurements, while the black squares 
are the best estimates after the correction (see the text). The size 
of spatial binning (0".6) is shown with the bar on the right side. 
The dotted lines present the values at Ad = 0. The positions of 1 
and 2 are indicated with the dashed lines. The right-hand vertical 
axis ticks the corresponding velocity. The FWHM-line width is 
corrected for the instrumental broadening with the assumption of 
a Gaussian profile. 

position 2). These two spectra show an asymmetric hne 
with S and S^ which are comparable with the average 
values of z ~ 6.5 LAE (see ^2^ . Thus, both of these 
components reside at high-z. The bottom panel of Figure 
[7] presents line-center offset and line width as a function 
of distance along the DEIMOS slit. We measure the line 
center and width by Gaussian-profile fitting. Note that 
the NB921 (or Lya) surface-brightness distribution is 
not homogeneous within the slit. In fact, Figure [T] implies 
that, from position 1 to 2, the profile center would shift 
towards the direction of red spectrum (magenta arrows 
in Figure[T|). We estimate the biases raised by this profile 
inhomogeneity with the NB921 image which has a seeing 
size comparable to our DEIMOS spectrum. We measure 
changes of the profile's center and standard deviation 
within the slit as a function of slit position in the NB921 
image. We calculate the correction factors in wavelength 
based on these spatial changes of A^i?921 profile, and ap- 
ply these correction factors to the original measurements 
of line-center offset and line width. The bottom panel 
of Figure [7] shows the corrected values (filled squares), 
together with the original measurements (open squares). 
We confirm that those biases are not large enough to 
alter general trends in line-center offset and line width. 



Figure [7] indicates that the Lya-line center shifts grad- 
ually by ~ 30 km s^^ between these two components. 
Thus, these components are not sources of an identical 
object produced by lensing because of this line shift. It 
indicates that it is unlikely that they are gravitationally 
lensed objects. Moreover, following t he method intro- 
duced in §2.4.2 of iPindor et all (|2006f) . we estimate the 
minimum brightness of lensing galaxy at z < 4 to be 
K = 23.2 based on our source redshift of 6.6 and image 
separation of I'M. Figure [T] presents no nearby bright 
sources with K < 23.2 to which can be ascribed a lens- 
ing object. There remains a special case where a very 
red foreground object with a 3.6/im-band detection and 
no optical-NIR counterparts magnifies an inhomogencous 
LAE with spatially different magnification factors. How- 
ever, there are little chance coincidences in the precise 
alignment of such an extremely-red lensing object. We 
thus conclude that our LAE is not likely to be a lensed 
source, but an intrinsically extended object. The size of 
the extended Lya nebula is > 17 proper kpc at z = 6.595 
which is estimated from the major axis of the isophotal 
area in the NB921 image (> 3".l; ^24|) . 

3.3. Kinematics of Resonance Lya Line 

Because Lya is a resonance line strongly scattered or 
absorbed by gas and dust, it is difficult to determine 
the gas kinematics by Lya observations alone. However, 
it is known that Lya line is a useful probe of gas in- 
flow/ outflow and cosmic reionization with detail ed mod- 
eling (j^^keneFaD [2Q03 [Dijkstrairn] [lOOli) . First, 
we obtain the line width of frwHM = 251 ± 21 km s~^ 
from the DEIMOS spectrum in the extraction aperture 
of 1".2 X 1".0 (i.e. 6.5 x 5.4 proper kpc^) around the Lya 
source center which corresponds to Ad — in Figure [71 
The line width is corrected for the instrumental broaden- 
ing with the assumption of a Gaussian profile. The line 
width can be twice as large as this value, if a blue half 
of Lya is completely absorbed by IGM with no effects 
of the Lya damping wing. From further inspection of 
the line-center offset and line width along our DEIMOS 
slit (Figure [7]), we find that the line-center velocity of 
Lya increases by Aw ~ 60 km s~^ from east to west in a 
range of ^ 2" {D = 10 proper kpc). Note that this small 
velocity offset is larger than the sizes of their error bars. 
^^ On the other hand, there are no significant changes of 
line width beyond the sizes of their error bars, although 
our spectrum implies an increase by ~ 50 km s~^ from 
east to west. 

3.4. Stellar Population and Mass 

We carry out x^ fitting of stellar synthesis models 
to the SED of this object based on total fluxes at the 
observed-frame of 0.9— 8.0^ m (Tablc[T|). We use our best 
estimates of the total fluxes with the associated la error 
for all of photometry points including those below the 
detection limits. Since our z'-band photometry is con- 
taminated by Lya emission and Gunn-Peterson trough, 
we estimate the emission-free continuum magnitude at 
9500A, Too. 95. We obtain too. 95 = 25.18 ± 0.73 based on 

^^ Although the instrumental spectral resolution is i? '^ 3600 
corresponding to tipwHM ~ 80 km s~^, the uncertainties of line 
centering by Gaussian fitting is as small as ~ 5 — 10 km s"-"^ (See 
error bars in Figure [71. 



Ouchi et al. 



NB921- and z'-band photomet r y by the method simi- 
lar to that of IShimasaku et aTl ( 2006f ). which takes ac- 
count the_contributions of the Lya Hne and IGM absorp- 
tion (|Madat]iri995[ ) with the response curves of NB921 
an d z' filters. We u se th e stellar synthesis models 
of iBruzual fc Chariot! ()2003l ) with dust attenuation of 
ICalzetti et alj (|2000D . Applying models of constant and 
exponentially-decaying star-formation histories with sets 
of metallicity in Z = 0.02 — I.OZq, we search for the 
best-fit model in a parameter space oi E{B — V) = — 1 
and age= 1 — 810 Myr, where the upper limit of stellar 
age is the cosmic age at z = 6.595. First, we assume 
the constant star-formation with a fixed metallicity of 
Z = O.O2Z0. Figure [H] presents the SED and the best- 
fit models. We find that the best-fit model has a stellar 
mass of M* = 3.5l2:6 x 1O^°M0 (i.e. 0.9-5.0 x 10^°Mq) 
with a reduced x^ oi 0.84. Because very weak photo- 
metric constraints are given in the rest-frame near UV 
{r-^ 0.2 — 0.3/im) critical to resolving the degeneracy be- 
tween extinction and age, we obtain no meaningful mea- 
surements on extinction and stellar age within the ~ la 
error. For examples, the sets of allowed parameters are 
(E[B-V], age[Myr])= (0.0,810), (0.3,200), (0.6,29), and 
(0.9, 3). On the other hand, our estimate of stellar mass 
(0.9 — 5.0 X IO^^Mq) has a moderate reliabihty due to the 
determinations of precise spectroscopic redshift and the 
rest-frame optical photometry on which stellar mass pri- 
marily depends. Because the sets of allowed extinction 
and age parameters cancel out the variance of mass-to- 
luminosity ratio (Af/L), stellar mass is obtained with a 
much smaller error than extinction and age (Papovich 
et al. 2001). It should be noted that this object is the 
most distant spectroscopically-confirmed galaxy whose 
stellar mass is constrained. Next, we change the metal- 
licity to Z = 0.2 — I.OZq and star-formation history to 
an exponential-decay time scale of r = 1 — 100 Myr. The 
stellar mass is estimated to be 1.7 — 3.2 x 10^'^ Mq and 
0.4 — 4.8 X IO^^Mq for the best-fit values and the la error 
ranges, respectively. We find that the general behaviors 
of the fitting is the same, and that the best-fit values and 
the la error ranges of stellar mass agree with those from 
the first assumptions (0.9 — 5.0 x IO^'^Mq) within a factor 

o f ~2. 

iMobasher et al.l ()2005D report a very massive galaxy 
at 2;phot = 6.5 with a stellar mass of 6 x IQ^^AIq based 
on a photometrically -selected galaxy. However, it has 
been pointed out that IMobasher et all 's object may be a 
lower redshift {z ~ 2) dust-enshroud ed starburst galaxy , 
based on a reanalysis of its photo-z (jDunlop et al.ll200"7f l 
and the detection of polycyclic aroin atic hydrocarbon 
emissi on features (ICharv et al.l [20071 ). On the other 
hand, iLabbe et al.l (|2006f ) study z-dropout galaxies with 
a photo-z of z ~ 7 in the Hubble Ultra Deep Field 
(HUDF) , and find that these galaxies have stellar masses 
of 0.1 — 1 X IO^^Mq. This mass range touches the 
low est-mass limi t of our object (0.9 — 5.0 x IO^^Mq). If 
the ILabbe et all 's photo-z objects are real high-z galax- 
ies, our object is likely to be a more stellar-massive 
galaxy th an those found in t he small area of H UDF. 
Similarly, [E gami et al.' (|2005D . ICharv et al] (|2005l ). and 
ISchaererfc PcUo (2005) have estimated stellar masses of 
gravitationally-lensed galaxies at z ^ 7 behind AbcU 
clusters to be ^ 10* — IO^Mq, which is ^ 1 — 2 orders of 
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Fig. 8. — Spectral energy distribution (SED) of the giant LAE, 
Himiko. The squares represent the total fluxes of this object taken 
from Table [l] The open symbols are data points that are not used 
for the SED fitting, since these bands are contaminated by the 
strong Lya emission line. The cyan, blue, magenta, and red lines 
show the best-fit stellar synthesis models with E{B — V) = 0.0, 
0.3, 0.6, 0.9, respectively, under the assumptions of constant star 
formation with a fixed metallicity of Z = 0.02^0. 

magnitude smaller than the stellar mass of our object. 

Because E{B — V) cannot be constrained, we can 
only obtain the lower-limit of star-formation rate of 
SFR > 34M0yr-i from the SED fitting, which is 
given in the case of E{B ~ V) = 0. This lower limit 
from the SED fitting is consistent with the SFRs esti- 
mated from the UV continuum of TO0.95 (SFR=25l;^2-^0 
yr~^) and from the Lya luminosity (SFR=36 ± 2Mq 
yr ~^) with no dust-ext inction corrections via formulae 
of iMadau. Pozzetti. fc Dickinson (1998) and Kennicutl 
(1998) + case B recombination, respectively. 



The SSFR of our object is SSFR > 1.6 x lO'^r 



9„,,-l 



at 



the stellar mass of 0.9 — 5.0 x 10^ Af©. In the plane of 
SSFR vs. stellar mass, our lower limit of SSFR is compa- 
rable to LBGs at z ~ 2 — 3 (see, e.g., Castro Ceron et al] 
l2008f ). However, the stellar mass of our object is, at 
least, one-order of magnitude larger than that of the av- 
eraged (stacked) LAEs at z ^ -^ 3 jNilsson et al.l 120071 : 
iGawiseret all 120071: iLai et al.l [20081 ). Our derived stel- 
lar mass is more similar to those of luminous LAEs 
at z = 3.1 — 5.7 that are bright enou gh to be iden- 
tified individually in infrared images (jLai et al.l l2007l : 
iFinkelstein et al.ll2008l : Ono et al. in prep ; cf. very faint 
LAEs in the HUDF bv lPirzkal et~al]|2007l ). 

We calculate four statistical measurements from our 
object's SFR, stellar mass, and number density listed in 
Table O and compare with those obtained by the other 
studies. The comparisons are useful to check how our 
object plays a role in the average volume of the Universe 
at z '^ 7. (i) We estimate the lower limit of UV luminos- 
ity function (LF) to be > 1.2 x 10"'^ mag'^ Mpc'^ at 
Afuv — -21.3. This limit is consistent with z ^ 7 UV 
LF of iBouwens et al.l (|2008f) . (ii) The lower Umit of cos- 
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mic SFR is > 4.3 x 10"'' M© yr"^ Mpc"^, which is one to 
two orders of maKnit ude smaller than that obtained by 
iBouwens et al.l (l2008fl. Thus, this lim it is consistent with 
the estimate of IBouwens et al.] (|2008[ ). The small contri- 
bution to the cosmic SFR leaves the possibility that sig- 
nificantly luminous objects like the one we discuss here 
could not be the major contributors of cosmic reioniza- 
t ion at z ~ 7 similar to at z '^ 6 (jYan fc Windhorst 
l2004f l. (iii) The lower limit of stellar-mass function is 
logp, > -5.9 Mpc-3 dex-i at log(M,) ~ 10.5. This 
limit is much low er than stellar mass function of z ^ 5 
dropouts given bv lMcLure et al.l ()2008l ). It indicates that 
our lower limit is compatible with the scenario of hier- 
archical structure formation, because a similarly massive 
system is more abundant at the recent epoch of z ~ 5 
than at z ~ 7. (iv) The lower limit of stellar-mass density 
is > 4.4 X 10"^ Mq Mpc"'^, which is consiste nt with that 
obtai ned in the HUDF (1.6xl0^M,7. Mpc'^^ lLabbe et all 
|2006[ ). All of these four statistical measurements corre- 
sponding to our object fit in the average properties of 
the Universe at z ~ 7. 

4. DISCUSSION 

4.1. Comparisons with Lya Blobs at z ^ 3 

We compare properties of our giant LAE with those of 
Lya blobs found at z ~ 3. Because objects experience 
more severe surface-brightness dimming at z = 6.6 than 
at z ~ 3, careful comparisons are needed. We p r oduce 
simulated 7V^9 21 images of blob 1 oflSteidel et all (|2000D 
and blob 28 of iMatsuda et all ()2004D redshifted to z = 
6.595 based on the narrow-band images of Matsu da et al.l 
11004). We choose the blob 1 and 28, because they 
have the brightest Lya luminosit y and the high est sur- 
face brightness, respectively, in IMatsuda et all 's cata- 
log. We carry out the simulations in the same man- 
ner as lSaito et al.l (j2006f ). but with an improved random 
noise whose amplitude exactly matches to those of real 
NB921 image. Figure |9| presents the simulated A^_B921 
i mages, together with t he original narrow-band images 
of IMatsuda et all (|2004f ). The simulated image of blob 
1 indicates that our observations would miss ~ 95% of 
blob I's L(Lya), and that no diffuse Lya nebula could be 
identified. The simulated blob 1 has (iVB921 [total], Ai^o, 
{SB)) = ( 25.5 mag, 0.8 arcsec^, 26.7 mag arcsec"^). 
We plot this simulated object in Figure JH The simu- 
lated blob 1 is indistinguishable from the cloud of nor- 
mal LAEs in the planes of ^iso vs. iV_B921 and Ai^o vs. 
(SB). On the other hand, the simulated blob 28 shows 
(A^B921 [total], Also, {SB)) = ( 24.1 mag, 3.7 arcsec^, 
26.3 mag arcsec"^). This simulated object is recogniz- 
able in Figure m However, this simulated object is more 
similar to the other LAEs with the second and third 
largest Also than to our object. There exist remarkable 
differences in Ai^o and {SB) (or NB921 magnitude) be- 
tween our object and the simulated objects. In other 
words, the large Aiso and the high {SB) of our object 
at z = 6.6 cannot be realized even for blob 28 whose 
very high surface brightness would minimize the effect of 
cosmological surface brightness dimming. It implies that 
our object might be a population that has not yet been 
identified at z ~ 3. 

Although the Lya brightness and profile of our object 
seem different from those of Lya blobs at z ~ 3, the other 
properties of our object show both similarities and differ- 
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Fig. 9. — Original narro w-band images of blob 1 and 28 at 
2 = 3.1 from iMatsuda et al . (2004) (top panels) and simulated 
narrow-band images at z = 6.6 (bottom panels). The size of each 
panel is 25" X 25". The intensity contours are presented with 
yellow lines. The contours represent (2, 3, 5, 10, 15) a levels of sky 
fluctuations, where the 2cr sky fluctuatio ns correspond to 28.0 and 
26.8 mag arcsec"^ in the original images IMatsuda et aH2004D and 
the simulated images, respectively. Note that the absolute values 
of contour levels in the simulated images are the same as the yellow 
contours in the NB921 image of Figure [T] In the original images, 
blob 1 is the large diffuse extended source, while blob 28 is the 
compact source at the center. Note that the original and simulated 
images include foreground and background sources neighboring the 
blobs in the 25" X 25" areas. 

ences. The line width of our object is comparable to some 
of the extended LAEs presented in Matsuda et al. (200q) 
in the plane of FWHM(Lya) vs. isophotal area, although 
the threshold of isophotal area is different from ours. On 
the other hand, there are no objects with such a lumi- 
nous Lya line in the range of ufwhm = 100 — 900 km s~^ 
in the catalog of Saito et al. (200_S). The velocity width 
of our object is 25 1 km s~ ^ which is significantly smaller 
than that of .Steidel et al.l 's blo b 1 with a broad-velo city 
wid th, - 1000- 1500 km s'^ (jOhvama et al.ll2003t see 
also lBower et al.ll2004l) . The velocity width of our object 
is closer to the one of blob 28 (wfwhm — 362 km s"^; Y. 
Matsuda in preparat ion). The stel lar masses of the sub- 
components shaping ISteidel et all's blobs 1 a nd 2 range 



from 10^° to 10"M,7^ ()Uchimoto et aLll2008l) , whic h is 



comparable to ours. Similarly, ISm ith et al.l ()2008l ) re- 
port that the stellar mass of their Lya blob at z = 2.83 
is as massive as 3 — 4 x 10^^ Mq. These stellar masses are 
comparable to ours within an order of magnitude. The 
estimated number density of our object is only 1.2 x 10~^ 
comoving Mpc ~'^, which is low er than those of z ~ 3 Lya 
blobs found bv lMatsuda et al.l (|2004) (3 x lO"" Mpc'^), 
ISaito et al.l (|2006[ ) (1 x 10"^ Mpc^^), and Yang et al. 
in preparation (3 x 10~^ Mpc""^). This number den- 
sity is also lower than the upper limit of number den- 
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sity of ISteidel eralT s two Lya blobs (2 x 10"^ Mpc"^), 
which are estim ated from the number of objects (2) and 
iMatsuda et all s survey volume (1.3 x 10^ Mpc^). Al- 
though these number densities depend on the criteria 
of sample selections and observing fields, our object at 
z = 6.6 would be as rare as (or even rarer than) these 
Lya blobs at 2 ^ 3. Regarding the environment, our gi- 
ant LAE reside s in a h igh density region of LAEs f ii2.1[) . 
IMatsuda eraD (|2004f ) find that, at z = 3.1, the dis- 
tribution of their Lya blobs trace the dense regions of 
LAEs. The similar spatial correlation between our ex- 
tended LAE and compact LAEs is also seen at z = 6.6. 
We also compare our objec t with Lya nebul ae around 
high-z radio galaxies fHzRGs: lMilev fc De Bre uck 2008). 
The Lya luminosity of our object is c omparable to those 
of HzRGs (- 10'^3_io44 erg s-^: [McCarthv et al. 1991 ). 
On the other hand, our object has a much smaller 
FWHM line wid th of Lya, 251 ±21 km s-\ than HzRGs 
(> 1000 km s~^: lMcCarthv et al.lll996l ). However, recent 
studies have found an HzRG with a strikingly relaxed 
dyn amics whose Lya FWHM line width is < 300 km 
s-i (jVillar-Martm et all [20071) . It would be difficuh to 
distinguish between our object and such a dynamically- 
relaxed HzRG powering a Lya halo based on dynamics, 
although our object has a radio flux fainter than that of 
a typical radio-loud AGN (Section 13. ip . 



4.2. Inferred Kinematics 

If this object forms a single virialized system whose 
possible two z'-band components would be just bright 
Hii regions in a disk (see §4.3|) . the dynamical mass 
from the rotation is naively estimated to be M-^at sin(j) 2± 
1 X IO^Mq by M,.otsin(i) = ([u^ sin(i)]2r)/G with Vc ^ 
Aw/2 = 30 km s^^ and r = D/2 = 5 kpc {^M, where 
i and G are the inclination of a rotating disk and the 
gravitational constant. On the other hand, from the 1- 
dimcnsional velocity dispersion of a„ = wfwhm/2.35 = 
107 km s~^ in a half size of the extraction aperture, R ~ 
6/2 = 3 kpc f iJ3.3p . we obtain the mass of random mo- 
tion of Mrand = 4x IQIOMq via Mrand = (5 / 3) (Sa^,) R / G , 
assuming a uniform sphere. If the blue half of Lya is 
absorbed by the external IGM, the velocity dispersion 
and mass are ~ 214 km s~^ and Mrand = 2 x 10^^ Mq, 
respectively. Thus, Afrand is 1-2 order(s) of magnitude 
larger than Af^ot in a reasonable range of inclination 
{i = 10 — 90°). If we consider the smaller radius for 
the Afrand estimate {R = 3 kpc) than that for the Afjot 
estimate (r = 5 kpc) , the difference of these dynamical 
masses defined in a common radius becomes even larger. 
If Lya reflects dynamics, our giant LAE would be a sys- 
tem more dominated by random motion than rotation. 
Note that these estimates of dynamical masses depend 
on the size of spectrum extraction aperture with the un- 
certainties of seeing smearing, and that these results are 
only true under the assumptions of the single-virialized 
system and no significant effects of resonant scattering. 

If this giant LAE is an outflow object whose Lya emis- 
sion is produced by shock heating (cf. iDiikstra fc Loebl 
|2008[ ) , the dynamical time scale required to form the ex- 
tended Lya nebula is ~ 7 x 10^ yr, where we assume the 
size of the major axis (~ 17 kpc) and the typical velocity 
width (upwHM — 251 km s^^). Since this time scale is 
as long as the cosmic time between z = 6.595 and 7.1, 
it would start making an ionized-bubble since z = 7.1 in 



this shock heating scenario. 

4.3. Nature of the giant LAE 

The nature of the Lya nebula of our object is not yet 
clearly understood within the currently available obser- 
vational data. There are five possible explanations: (1) 
halo gas photoionized by a hidden AGN, (2) clouds of 
Hii regions in a single virialized galaxy, (3) cooling gas 
accreting onto a massive dark halo associated with an 
initial onset of starburst at the halo center, (4) merging 
bright LAEs with clouds of Hii regions, and (5) outfiow- 
ing gas excited by shocks or UV radiation from starbursts 
and/or mergers. There is a chance of (1), although no 
positive evidence of an AGN is found. Our object shows 
the lack of Nv line, no line broadening, and no detections 
in X-ray, MIR, submm, and radio bands (Figure[6|). But, 
these constraints are not strong enough to discard the 
possibility of a hidden AGN. The case of (2) seems sur- 
prising, because it means that such a large galaxy exists 
in a very early epoch of z = 6.6. This galaxy would have 
a size of > 17 kpc and possibly two large star- forming re- 
gions in a disk seen in our z' band image ( §3.2|) . However, 
there is a chance to explain this large galaxy in the frame- 
work of Cold Dark Matter (CDM) models. We estimate 
properties of the most massive dark halo whose number 
density is the same as o ur object at z = 6. 6 . Bas ed on the 
analytic CDM model of lSheth fc TormenI ()1999l ). we find 
that the dark halo has a radius of 47 kpc and a mass of 
1 X 10^'^ Mq with a circular velocity of 380 km s~^ All of 
these values of the dark halo are significantly larger than 
those of our giant LAE measured via radiation, indicat- 
ing that the single-galaxy picture can be compatible with 
the CDM model in terms of halo properties. But it is 
not obvious that such a big virialized baryonic system at 
z = 6.6 can be reproduced in the scheme of CDM model. 
The case of (3) is possible with the potentially large EWq 
of 57 — 402A. Moreover, our DEIMOS spectrum presents 
a possible weak red peak at ^ 9245A (Figure O which 
is redshifted from the main Lya peak (9233A) by ~ 400 
km s~^. This possible red peak may be similar to the 
one produced by IG M infalling into a colla psing cloud, 
which is claimed by IDiikstra et al.l (|2006a !,rp) . On the 
other hand, the comparable SFR values from Lya and 
UV continuum (i i3.4p may not prefer this scenario, since 
the Lya luminosity could be explained solely by normal 
star- forming activities that the UV continuum indicates. 
However, this argument is not strong due to the under- 
estimation of SFR from Lya luminosity. In fact, the 
intrinsic Lya luminosity could be brighter by a factor of 
20 if our giant LAE has diffuse Lya components simi- 
lar to ISteidel et al.l 's blob 1 ( §4.1[1 . The explanation of 
(4) would be reasonable, given the possible existence of 
two peaks in z' band with a separation of l".l (6.0 kpc). 
The A^i?921-image profile can be reproduced by two ex- 
ponential disks with a reasonably small half-light radius 
of Rhi = 0".3 corresp onding to Rhi — 1.6 kpc f ij3.21 
cf. ISimard et al.lfl999[ l. The merger would induce star- 
formation activities, and could produce the bright Lya- 
line and UV-continuum emission. Figure |7| shows that 
velocity widths at position 1 and 2 are similar, implying 
that dynamical masses of these components would be 
comparable. We may be witnessing the site of a major 
merger near the reionization epoch. On the other hand. 
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our object shows a Lya nebula potentially larger than 
the isophotal scale with the major axis of 3".l (17 kpc; 
tj2.ip . If the Lya nebula really extends beyond the 2(7- 
level isophotal area, it becomes difficult to explain the 
Lya morphology with a profile of two merging LAEs. 
The (5) case seems plausible, since our object has a rel- 
atively high SFR and a large stellar-mass with possible 
multiple components in z' band fi )3.2p . In either case 
of (4) or (5), our object would be a massive galaxy in 
formation with significant star-forma tion contributing to 
cosmic reionization (cf. Illiev et aLl i2006 ) and/or with 
outflo ws for the metal enrichment of IGM (jBouche et al.1 
|2007| ). Since our object has the small velocity offset 
{Av = 60km s~^) and the line width (wfwhm = 251 
km s~^), the dynamics of merger or outflow would have 
to be well collimated to the direction perpendicular to 
the line of sight. 

The angular size of the Lya nebula is > 17 proper kpc 
( §3.2p . which is comparable to the diameter of the stellar 
disk of the present-day Milky Way. It is impressive, if we 
consider that the age of the Universe at z = 6.595 is only 
6% of the one of the present-day Universe. Moreover, 
such an extended Lya source is very rare in the cosmo- 
logical volume only with the number density of 1.2 x 10~^ 
comoving Mpc~'^ at z = 6.6 (' i)2.2p . If our selection of 
large Lya nebula does not miss a signiflcant fraction of 
massive galaxies at this early epoch (2 — 6.6), our object 
could be an ancestor of a bright-cluster or cD galaxy, and 
should be a good laboratory of massive-galaxy formation 
near the reionization epoch. 

4.4. Future Prospects 

In section 4.3, we have found that the currently avail- 
able data do not provide a clear answer to the question 
about the nature of this object. It is obvious that deeper 
NIR and infrared images of this object can be taken with 
Hubble and Spitzer Space Telescopes to constrain SFR, 
dust extinction, and stellar age, which will trace back 
through the star- formation histor y of this object (e.g. 
lYan et al.l [200I lEvles et al.ll2007[ ). These deep images 
are critically important to give stronger constraints on 
the SED of our object (Figure [5]) that characterizes SFR 
and dust extinction, and to identify Lya photons that 
are not originated from star-formation activities but from 
others, such as cold accretion. The high-resolution NIR 
image of Hubble Space Telescope is useful to investigate 
the possibilities of mergers and outflows. A J-band spec- 
trum will test the existence of AGN with a relatively 
strong Civ emission line. We estimate an expected Civ 
flux to be '^ 2 and 4 x 10~^^ erg s~^ cm~^ for Seyfert II 
and QSO, assuming the ratios of few / fhya ^ 0.2, and 
fcw/fhya ~ 0.5, respectively (lMcCarthvlll993f ). These 
flux limits are achieved by deep NIR spectroscopy with 
8m-class telescopes. Moreover, we can characterize star- 
formation activities and metal enrichment in our Lya 
nebula by deep submillimeter and millimeter observa- 
tions with Atacama Large Millimeter/Submillimcter Ar- 
ray (ALMA) which will start the operation in 2012 pre- 
ceded by the early-science operation. ALMA observa- 
tions will allow us to investigate emission from dust and 
molecular-clouds in our object. A detection of spatially- 
extended metal line from the Lya nebula could reject the 
possibility of cooling accretion of primordial gas, and pro- 
vide an independent probe of star-formation properties. 



We estimate an expected intensity of dust emission in 
850^m, S'(850^m), and a flux of molecular CO(6-5) fine, 
S{CO), assuming typical parameters of local starbursts. 
We start the calculations from our lower-limit of SFR, 
34M0yr~^. The far-infrared luminosity, L(FIR), is cal- 
culated from SF RlMc^yr-^ = 1.7 x 10-^° L {FIR) [Lq] 
(|Kennicuttlll99l . We obtain L(FIR) = 2.0 x lO^L©. 
Assuming the modified blackbody radiation with a dust 
emissivity index of /? = 1.3 arid a d ust temperature of 
Tdust = 35.6 K (jDunne et al.l l2000f) . we estimate the 
850^m dust emissi on from L(F IR) to be 5(850/im) — 
0.28 mJy (see, e.g.. lOuchi et al.l lT999). This moderately 
bright 850^m emis sion is expected b ecause of the neg- 
ative k-correction (jBlain et al.l l2002f ) . Again from the 
-L(FIR) value with the same modified blackbody radi- 
ation, the dust mas s is Mnnst, = 6.0 xJ^Mq via the 
relation presented in iDe Breuck et al.l (|2003[ ) . The mass 
of molecular hydrogen is M{H2) = 3.0 x lO^M© which 
is calculate d with the relation of M(H2)/Mdust = 50 
f^ 25 - 75: ISeaauist et al.l l200l:'^ ^g et al.lllBl . Fi- 
nally, the flux of molecular CO (6-5) line is S'(CO) = 0.1 
Jy km s~^ with the assumptions of an H2-to-C O con- 
version factor of 0.8 (jDownes fc Solomonlll998D. and a 
line r atio of C0(6 - 5)/C0(l - 0) = 0.5 (iBavet et~all 
120061 ). In summary, we expect S'(850/im) > 0.28 mJy 
and S{CO) ^ 0.1 Jy km s~^, considering that our SFR 
is the lower hmit (> 34M0yr~^). We use ALMA Sensi- 
tivity Calculator ^', and estimate the on-source integra- 
tion time to be < 0.2 and < 3 hours for ba detections 
of a 850/im continuum and a CO (6-5) line, respectively. 
Here we assume the large beam size of 1" in a compact 
configuration of 50 12m-arrays for the point-source detec- 
tion, the band width of 16 GHz, and the CO-line width 
of 250 km s-i (jNishivama fc Nakail[200l[ l with a 50 km 
s~^ spectral resolution. Either of an 850/im-thermal con- 
tinuum or a molecular-CO(6-5) line may be detected in 
reasonable observing time under the assumptions of local 
starbursts. If our object does not have dust or molecular 
gas as much as the local starbursts, a deficit of dust or 
molecular-line emission would be identified by ALMA ob- 
servations. In either case, dust and molecular-gas prop- 
erties of our object could be characterized in a few years. 
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